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Abstract 



'sj" ■ To address the issue of whether tri-bimaximal mixing (TBM) is a softly-broken hid- 

^^ . den or an accidental symmetry, we adopt a model-independent analysis in which we 

f~^ \ perturb a neutrino mass matrix leading to TBM in the most general way but leave 

the three texture zeros of the diagonal charged lepton mass matrix unperturbed. 
We compare predictions for the perturbed neutrino TBM parameters with those ob- 
tained from typical 50(10) grand unified theories with a variety of flavor symmetries. 
C^ ' Whereas 50(10) GUTs almost always predict a normal mass hierarchy for the light 

neutrinos, TBM has a priori no preference for neutrino masses. We find, in particular 
for the latter, that the value of \Ue2,\ is very sensitive to the neutrino mass scale and 
ordering. Observation of |C/^s3p > 0.001 to 0.01 within the next few years would be 
incompatible with softly-broken TBM and a normal mass hierarchy and would sug- 
gest that the apparent TBM symmetry is an accidental symmetry instead. No such 
conclusions can be drawn for the inverted and quasi-degenerate hierarchy spectra. 
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Neutrino oscillations seem to point towards a lepton flavor structure completely dif- 
ferent from that of the quark sector. In particular, the PMNS lepton mixing matrix has 
a very different structure from that of the CKM quark mixing matrix. Nevertheless, this 
seemingly incompatible feature can be reconciled in Grand Unified Theories (GUTs), where 
quarks and leptons belong to the same multiplets. In particular, GUTs based on SO{10), 
which allow for a seesaw mechanism [1] without adding singlets by hand, have been fre- 
quently studied in the past ten years; see e.g. [2]. The models often specify in addition a 
particular fiavor symmetry with charges assigned to the fermion and Higgs 5*0(10) mul- 
tiplets, although so-called "minimal" Higgs models may rely on no fiavor symmetry at 
all. 

Another more recent approach for explaining the neutrino mixing scheme has involved 
the introduction of lepton fiavor symmetries. The goal of such models has been to reproduce 
the approximate tri-bimaximal mixing (TBM) form observed by Harrison, Perkins and 
Scott among others [3]: 
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Here P = diag(l, e*", e*^) contains the two Majorana phases, a and (3, whereas the Dirac 
phase S remains unspecified. Lepton fiavor symmetries such as A4, 6*4, and S3 have been 
applied which lead to the TBM in a straightforward manner [4]. Even more recently 
symmetries such as T' have been introduced as an extended fiavor symmetry in order to 
treat the quark sector as well in a self-consistent fashion [5]. 

The question then arises whether the TBM symmetry is a presumably softly-broken 
hidden symmetry, or whether it is an accidental symmetry in nature. Note that the fiavor 
symmetries originally introduced with 5*0(10) models were not designed to reproduce 
the TBM matrix per se, but rather were designed to reproduce quark and lepton mixing 
schemes in approximate agreement with the then known mixing data. Even with more 
refined data and fits to the data now available in the literature, many 5*0(10) models have 
still survived. For reference, we quote the current best-fit values and la (3o") ranges for 
the neutrino oscillation parameters as given in [6]: 
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For exact TBM, the mixing angles correspond to 



sin^^ 



12 



1 

3 



sin^ e 



23 



1 
2 



sin^^ 



13 



. 



(3) 



In order to address this issue and provide a partial answer to the question raised, we 
shall study (what we consider) reasonable deviations from TBM and compare those results 
with the predictions of S'O(IO) models available in the literature. In doing so, we will not 
assume a particular model leading to TBM, but instead take the corresponding neutrino 
mass matrix for exact TBM at face value after adopting a basis in which the charged 
leptons are real and diagonal, i.e., Ui = 1. Let us stress here the following: we take the 
point of view that some unknown flavor symmetry generates TBM and work in the charged 
lepton basis. The flavor symmetry results obtained in any other basis can readily be cast 
into this form, so our results will apply in general. With the specified choice of basis the 
mass matrix, m^, uniquely giving rise to TBM is 
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Here ??i^'^^ = diag(mi, m2, ms) and the parameters A,B,D are in general complex and 
functions of the neutrino masses and Majorana phases: 
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In this short note we will investigate in a general manner deviations from the tri-bimaximal 
mixing texture. Our Ansatz is to modify the structure of the mass matrix by multiplying 
each element of Eq. (jl]) with an individual complex correction factor, e,: 
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Here the complex perturbation parameters are taken to be |ej| < 0.2 for z = 1 — 6 with 
their phases (pi allowed to lie between zero and 27r. 

Note that had we chosen instead to perturb the original three parameters A 
D with complex parameters e^i, e^ and eo-, the neutrino masses, m 
altered but the mixing matrix would remain TBM. Instead we perturb the neutrino mass 
matrix as above but demand that the three texture zeros in the diagonal charged lepton 
mass matrix, m^, remain unperturbed. The same perturbation prescription then applied 
to 7711 simply results in a diagonal phase transformation acting on the matrix [/^, which 
can be rotated away in t/pMNS = U} Utqm- 



Still one may insist on applying corrections from the charged lepton sector. For example, 
with [/pMNS = Uj U,y, one can assume that U,^ corresponds to tri-bimaximal mixing and 
that the correction is given by 



f/, ~ -A 1 . (7) 




Then in the basis in which the charged leptons are diagonal, the neutrino mass matrix 
reads 

ml = Ul {my)TBM Ul . (8) 



In the case of a normal hierarchy where /I ~ i? ^ JAniQ/S ^ Z^ ~ JAm\, one finds 
for the e/i entry that {m'^)e^ ~ 5 (1 — ^ |). Since § | — | jArnjJAm^X < 9 A, we 



require A ^ 0.02 in order to have a soft-breaking perturbation to (mj,)TBM of less than 



20%. A similar number holds in the inverted hierarchy, where for a = tt/2, A :i^ J Am\/3 



and B ^ —2J Am\/3, the yU/i entry receives the largest corrections, roughly 8A. Again, 
demanding a soft-breaking correction less than 20% requires A ^ 0.025. It is then easy to 
see from [7] with f/pMNS = ^/ ^tbm and Ui from Eq. (jTj), the following small deviations 
from TBM are obtained: 

\Ue3\' = y and sin^ 2^23 ^l-^\' . (9) 

The implied small values of A lead to [t/^ssP well below 10~^ and sin^2^23 very close to 1. 
In the spirit of this note we do not tolerate a leptonic correction with as large a value as 
A = 0.22 for the soft breaking perturbation. We prefer to work in the charged lepton basis 
and proceed as indicated in the previous paragraph. 

Radiative corrections also lead to perturbations of a tri-bimaximal mass matrix [8]. The 
relevant small parameters depend on the charged lepton masses, so that the r contribution 
is enough to consider. In this limit the er and fir elements are multiplied with (1 + er), 
while the tt entry is multiplied with (1 + 2er)- The small parameter is defined as e^ = 
c Ye^ In ^, where v = 174 GeV and c is given by 3/2 in the SM and by -(1 + tan^ /5) 
in the MSSM. Demanding that |2e^| be less than 0.2 leads for Mx = 10^^ (10^) GeV to 
tan/3 ^ 71.4 (97.1). Therefore, our corrections include radiative corrections up to this 
huge value of tan f3. A detailed analysis of radiative corrections to TBM can be found in 
Ref. [9,10]. 

Returning to the perturbed matrix in Eq. (Q, we will vary the complex e parameters, 
diagonalize the resulting m,^'s and study the predictions for the neutrino mixing angles. 
The results obtained for the perturbed mixing matrix will of course depend on the neutrino 
mass ordering and scale. We shall find that \Ue3\ depends most sensitively on these observ- 
ables. Luckily, one expects a sizable improvement on its current upper limit of If/esI ~ 0.2 
in the near future. We then compare our findings from broken TBM with predictions 



from successful SO (10) GUTs and study how one may distinguish these two approaches 
experimentally. 

Let us first consider normal mass ordering. The strategy we adopt is as follows: we 
fix m^ = 0.050 eV, and take as starting values for the other masses, m2 = 0.0095 eV and 
mi = 0.0037 eV, corresponding to Aml^ = 7.66 x 10"^ eV^ and Amj^ = 2.49 x 10"^ eV^ 
values well within the center of the currently allowed region. We shall allow for a 20% 
variation around the initial values of m2 and mi and vary the phases a and (3 between 
zero and 27r. Furthermore, the complex perturbation parameters in Eq. ([6]) are also varied 
within |ej| < 0.2 for each i = 1 — 6, with the full range of phases allowed for each. For 
each choice of parameters the resulting mass matrix is diagonalized and, if the outcome 
is within the current 3a range from Eq. ([2]), the point is kept. Note that with a maximal 
perturbation of 20% of the individual mass matrix elements, two entries can have a relative 
variation of 40%, which is quite generous. 

The scatter plots in Figs. [T] and [2] show the results of this analysis. We see that |f/e3p 
is predicted to be rather small and lie below 10"'^. This number should be compared 
with the expected sensitivity of the Double Chooz reactor experiment [11], which will 
start data taking in 2009, and will reach a 90% CL. limit of 0.018 after one year with 
one detector, and 0.005 after 3 years of operation with both detectors. The Daya Bay 
experiment [12], presumably starting after Double Chooz, is expected to improve the limit 
by another factor of two. Our results show that neither of the two experiments is expected 
to find a positive signal, if a flavor symmetry predicting tri-bimaximal mixing and a normal 
mass hierarchy is broken by less than 20%. The same is true of course for the currently 
running long-baseline MINOS [13] and OPERA [14] experiments. It will take the first 
generation superbeam experiments, or perhaps even more advanced technologies such as 
/3-beams or neutrino factories, to probe |t/e3p in the range below sin^ ^13 = 10~^ predicted 
by the perturbed TBM results. On the other hand, sin^ 612 is uniformly populated over its 
experimentally allowed range and appears to have no bearing on the issue raised. 

For comparison, we also give in Tabled] the results for thirteen 5*0(10) GUT models, 
all of which involve a conventional type I seesaw mechanism and predict a normal mass 
hierarchy for the light neutrinos. All of these models predict all three mixing angles in 
their currently allowed range. According to the names of the authors, the references are 
A [15], AB [16], BB [17] BM [18], BO [19], CM [20], CY [21], DMM [22], DR [23], GK [24], 
JLM [25], VR [26] and YW [27]. For details, we refer to the cited works, and also to recent 
model compilations [28-30]. Note that the S'O(IO) predictions are well separated from 
the results of perturbed TBM and are accessible or more nearly accessible to the reactor 
experiments discussed above. The only exception is the GK model. We note in this respect 
that this model is very much challenged by its relatively large predictions for lepton flavor 
violating decays like fi ^ e-y [29]. 

It is of interest to present some approximate analytical results to support the numerical 
work leading to the scatter plot in Fig. [H We will not try here to diagonalize the fully 
perturbed mass matrix, but rather estimate the implied order of magnitude of jf/esl and 
sin^ 2623. A general statement, independent of the mass ordering, is that if the perturbation 
occurs only in the ee or fir entry of the tri-bimaximal rrii, from Eq. (jlj), then the resulting 
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0.0024 


0.30 


GK [24] 


NH 


1.00 


0.00059 


0.31 


JLM [25] 


NH 


1.0 


0.0189 
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Table 1: 5*0(10) models and their predictions for the lepton mixing angles. If ranges are 
given we take the central value. Also given are the constraints, if any, on the mixing angles 
for the three possible mass orderings from the softly-broken tri-bimaximal mixing mass 
matrices. 



mass matrix is still ^-t symmetric. Consequently the relation \Ue3\ = 1 — sin^ 2^23 = 
still holds in this case. The largest deviation of {Uesl from zero occurs when the e/i and 
er entry of the tri-bimaximal rrii, are perturbed [10,31]. If we set all other perturbations 
to zero, the extreme case occurs when the e/i element is multiplied by (1 — e) while the er 
element is multiplied by (1 -|- e), where e is here real. We can then diagonalize Eq. ([6]) and 
find (ignoring the phases and neglecting rui for simplicity) that 
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where R = Amg/Am^ is the ratio of the solar and atmospheric mass-squared differences 
with ArriQ = ml —ml and Am,\ = \m,l — ml\. This is actually quite close to the numerical 
result, and the small discrepancy can be explained by the effects of non-zero mi and the 
other small terms including e^ . The smallest value for sin^ 26'23 is achieved for a pertur- 
bation in the /i-r block of the tri-bimaximal m,^, [10,31]. The extreme case occurs for a 
multiplication of the /x/i entry by (1 + e) and of the tt entry by (1 — e), and yields 



sin2 2^„ ~ l-e^ >0.96 . 



(11) 



This is also quite close to the numerical result. 



For the inverted hierarchy case, we note that such a stable hierarchy is difficuh to 
obtain in SO{10) models which do not have a type II seesaw structure, i.e., if there is 
no direct left-handed Majorana contribution arising from a real or effective Higgs triplet. 
The 5*0 (10) model from Ref. [32], which has a negligible type II (triplet) contribution, 
is able to fit an inverted hierarchy, but is not very predictive in what regards the mixing 
angles. The model from Ref. [33] concentrates on embedding a neutrino mass matrix with 
vanishing diagonal elements, assumes type II (triplet) dominance, and has a best-fit of 
|f^e3p = 0.0025 in case of an inverted hierarchy but no other mixing angle predictions. 

The results for perturbed tri-bimaximal matrices arising from an inverted hierarchy are 
shown as scatter plots in Figs. [3] and HI Here we have started with fixed m2 = 0.05076 
eV, mi = 0.050 eV and m^ = 0.0114 eV, corresponding to Am'^i = 7.66 x 10~^ eV^ and 
Arrig^ = —2.37 x 10~^ eV^. Proceeding as in the normal hierarchy case we have varied 
the phases, a and (3, and masses mi and m^ within 20% of their starting values. To guide 
the eye, we have drawn with a dashed line the value of |f/e3p = 0.001, which is roughly 
the upper value found in case of a normal mass hierarchy which separates the perturbed 
results from 5*0(10) GUT predictions. It is clear from the plots in Figs. 3 and 4 that |f/e3p 
can easily be around, and even above, 0.01 and is therefore testable in up-coming reactor 
experiments, unlike the normal hierarchy case. Again sin^ 9i2 is unconstrained. 

In Figs. [Hand Owe show results for restricted values of a = and it/ 2 appearing in the 
parameters A and B of the TBM neutrino mass matrix. We find that the largest values 
of |f/e3p typically occur if the phase a is around tt/2. This value implies that the effective 
mass governing neutrinoless double beta decay (Oi//3/3) takes its minimally allowed value: 



m) :i^ JAm,\ cos 2^12 — | JAm,\. For the normal hierarchy case, (m) is expected to be 



smaller still. On the other hand, if the phase a is zero or vr, (m) ~ JAm,\, and |t/e3p is 
tiny. We also note that 623 can deviate more sizably from maximal mixing, if neutrinos are 
inversely ordered, and that the largest deviation occurs for a = tt/2 as shown in Fig. [61 



Turning to analytic estimates for inverted ordering, for a = 7r/2 one has A ~ \/Am,\/3 




and B ~ —2JAm,\/3. Consider first a perturbation of the e/x entry with (1 + e) and of 
the er entry with (1 — e) for real e. In this case 

< 10^2 and sin^ 2^23 ^ 1 - (^) e^ ^ 0.87 . (12) 

In the case of a = 0, we have B/A ~ | Am^/Am^ and find that |f/e3p is at most of order 
{eB/Ay ~ 10^^ and therefore completely negligible. In the case of perturbed /i/i and tt 
entries, sin^ 2^23 — 1 — e^ ^ 0.96. The agreement with the figures is quite reasonable. 

For completeness we study also the case of quasi-degenerate neutrinos. For definiteness, 
we consider the case of normally ordered neutrinos, choosing fixed ma = 0.1 eV, 7712 = 
0.08778 eV and mi = 0.08735 eV, corresponding to Amji = 7.59 x 10"^ eV^ and Amji = 
2.37 X 10"'^ eV^. The procedure is the same as before and the results are shown in Figs. [7] 
and [HI Plots for an inverted ordering of quasi-degenerate neutrinos look basically identical. 
As expected, both sin^ 2^23 and lUesl"^ deviate more sizably than before from their initial. 



tri-bimaximal values. Since this quasi-degenerate case is more similar to the inverted one, 
it is not surprising that the interplay of the Majorana phase a, the effective mass and 
the deviation from maximal ^23 and zero ^13 are similar. Analytically, one finds that an 
enhancement of roughly ml/Amj^ ~ 4 occurs for the upper limits of |f/e3p- All mixing 
angles are populated in their allowed ranges. 

In summary, we have raised the question whether the approximately tri-bimaximal 
mixing observed in the lepton sector results from a hidden symmetry or whether it is 
accidental. Early proposed mass matrix models based on 5*0(10) family symmetry were not 
designed to lead to tri-bimaximal mixing, although a number of them are still successful. 
To study this issue we have adopted a model-independent approach and perturbed the 
TBM neutrino mass matrix elements about their central values by 20% in the lepton flavor 
basis. The charged lepton mass matrix is trivially perturbed by keeping the off-diagonal 
three texture zeros intact. We found that |f/e3p and the neutrino mass scale and ordering 
are of importance in the problem. In general the value of sin^ 9i2 is not constrained, and 
precision measurements of its value will not help in settling the issue raised. The other 
mixing parameter, sin^ 2^23 has only limited impact. A most striking result obtained is 
that for a normal neutrino mass hierarchy, the predicted perturbed TBM values of |f/e3p 
lie below 10~^, while the type I seesaw SO{10) models typically predict values above this. 
For an inverted or quasi-degenerate neutrino mass hierarchy, on the other hand, only two 
of the studied 5*0 (10) models apply, while the allowed perturbed values of |f/e3p can range 
noticeably higher, up to the present experimental limit. An interesting correlation with 
the value of the effective mass governing neutrinoless double beta decay is observed. While 
the question posed is unanswerable at this time, we can conclude that observation of {Uesl 
within the next few years would be incompatible with TBM and a normal mass hierarchy. 
Clearly it will be necessary to determine both \Ue3\ and the mass hierarchy in order to 
address the posed question. 
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Figure 1: Scatter plot of sin^ 2^23 against |t/e3p for perturbed tri-bimaximal mixing and a 
normal mass hierarchy. Also given are predictions of thirteen 5*0(10) GUT models. 
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Figure 2: Scatter plot of sin^ ^12 against |f/e3p for perturbed tri-bimaximal mixing and a 
normal mass hierarchy. Also given are predictions of thirteen 50(10) GUT models. 
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Figure 3: Scatter plot of sin^ 2^23 against |t/e3p for perturbed tri-bimaximal mixing and 
an inverted mass hierarchy. 
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Figure 4: Scatter pfot of sin^ 9i2 against |f/e3p for perturbed tri-bimaximal mixing and an 
inverted mass hierarchy. 
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Figure 5: Scatter plot of the effective mass against |f/e3p for perturbed tri-bimaxiinal 
mixing and an inverted mass liierarcliy with two different choices of the Majorana phase 
a, with the a = 7r/2 cluster on the left and the a = cluster on the right. 
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Figure 6: Scatter plot of sin^ 26'23 against |f/e3p for perturbed tri-bimaximal mixing and 
an inverted mass hierarchy with two different choices of the Majorana phase a, with the 
upper cluster referring to a = 7r/2 and the lower to a = 0. 
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Figure 7: Scatter plot of sin^ 2^23 against |t/e3p for perturbed tri-bimaximal mixing and 
quasi-degenerate neutrinos. 
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Figure 8: Scatter pfot of sin^^i2 against |f/e3p for perturbed tri-bimaximal mixing and 
quasi-degenerate neutrinos. 



